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Amendments to Claims: 

This listing of claims will replace all prior revisions, and listings, of claims in the 
application: 

Listing of Claims: 

1 . (Currently Amended) A method for solving a Model Predictive Control problem to 
control a plant including the steps of: 

a) forming a Large Sparse Matrix Equation, LSME, based upon the Model Predictive 
Control problem; 

b) repeatedly, at least once per controller update, sampling commands and feedback 
sensors to pose a determination of actuator commands in terms of a solution of a quadratic 
programming problem based upon model predictive control; 

c) forming a square root of H, Hr, where H is a submatrix of the LSME that varies only 
once per controller update and Hr is a block triangular matrix in which each block along a 
diagonal is triangular; 

d) forming a square root of a Large Sparse Matrix of the LSME, LSMroot, using Hr in 
each of a plurality of iterations of a quadratic programming solver within each controller update, 
and LSMroot is block triangular matrix in which each block along the diagonal is triangular; and 

e) completing the solution of the LSME based upon LSMroot in each iteratio n; and 

f) outputting the plurality of actuator commands based upon the solution of the quadratic 
programming problem from the quadratic programming solver in each controller update to a 
plurality of actuators to control operation of a plant. 

2. (Original) The method of claim 1 wherein Hr is calculated in a first of the plurality of 
iterations and is not recalculated for subsequent ones of the plurality of iterations. 

3. (Original) The method of claim 2 wherein in said step c), Hr is formed without first 
forming H. 
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4. (Original) The method of claim 3 wherein said step c) is performed with a sequence 
of QR factorizations. 



5. (Original) The method of claim 3 where H = Hr'*S*Hr and Hr = L W" where 
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6. (Original) The method of claim 5 where the submatrices F, G, and M are produced 
by the recursive sequence of QR factorizations: 
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for various future controller update time points in the prediction horizon. 
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7. (Original) The method of claim 1 wherein the LSMroot is formed in said step d) 
without first forming the LSM. 

8. (Original) The method of claim 7 wherein said step d) is performed by solving a 
normal or generalized QR factorization. 

9. (Original) The method of claim 7 wherein said step d) is performed by solving a 
generalized QR factorization customized for an interior point method. 

10. (Original) The method of Claim 9 where the generalized QR factorization is 



where Pin and S are inputs and triangular P and, possibly, U are outputs, and where 



1 1 . (Original) The method of claim 7 wherein said step d) is performed by solving a 
generalized QR factorization customized for an active set method. 

12. (Original) The method of Claim 1 1 where the generalized QR factorization is 

s V2p = S V2 V . pin ^ where pin = L -TjT ^ jjT s mU = s 

where Pin and S are inputs and P and Sr, and possibly, U are outputs and where 



13. (Original) The method of Claim 12 where the generalized QR factorization is 
accomplished via a standard Householder transform based QR algorithm that is 
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modified by replacing the sequence of Householder transforms U ; = 




which are normally used with a corresponding sequence of generalized 

transformations U ; = \ I - — 1 - U is unneeded but, if desired, is the 
t v T S.v) 

cumulative product of the Ui. 

14. (Original) The method of claim 13 where the result of generalized QR factorization, 
which is a triangular matrix consisting of rows containing either purely real numbers 
or purely imaginary numbers, is further factored into the product of a diagonal matrix, 
S l R 12 , which has cither 1 or the square root of -1 on the diagonal, and a purely real 
triangular matrix, and wherein the square of this diagonal matrix, Sr, will then have 
either 1 or -1 on the diagonal, and Sr and P are purely real numbers. 

15. (Original) The method of claim 1 1 further including the step of factoring the LSM to 
use only real numbers as: 



where S and SR are constant diagonal matrices with only 1 or -1 entries along the 
diagonal, and where LSMroot is a block triangular matrix where each diagonal submatrix is also 
triangular. 

16. (Previously Presented) The method of claim 1 wherein the LSME has the form: 
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where z is a vector containing states, controls and state equation adjoint 
variables for each time point, grouped by time; m includes adjoint variables for inequality 
constraints, grouped by time; f and K are vectors; H and J are banded matrices; and T is 
a diagonal matrix or zero. 

17. (Currently Amended) A method for controlling a plant s ystem including the steps 
of: 

a) Receiving receiving a plurality of sensor signals indicating current conditions of the 
system; 

b) Receiving receiving a plurality of commands; 

c) Determining determining a desired dynamic response of the system based upon the 
commands; 

d) in each of a plurality of time steps, formulating a problem of achieving a desired 
dynamic response for a window spanning multiple time steps as a solution to a quadratic 
program problem using methods of model predictive control and the sensor signals; 

e) solving the quadratic programming problem in each time step using an iterative 
algorithm in which a Large Sparse Matrix Equation, LSME, is formed based upon the model 
predictive control problem and a quadratic programming algorithm; 

f) in a first iteration of the iterative algorithm, forming a square root of H, Hr; 

g) in iterations including the first iteration and iterations subsequent to the first iteration 
of the iterative algorithm, forming a square root of a Large Sparse Matrix of the LSME, 
LSMroot, using Hr in a quadratic programming solver; and 

h) in iterations including the first iteration and iterations subsequent to the first iteration 
of the iterative algorithm, completing the solution of the LSME based upon LSMroot; 

i) generating a plurality of actuator commands based upon the solution to the quadratic 
programming problem from said step e); and 

j) outputting the plurality of actuator commands to a plurality of actuators acting upon the 
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18. (Currently Amended) A model predictive control system for controlling a plant 
comprising: 

a plurality of sensors indicating a current state of the system; 

a desired trajectory generator for creating a desired dynamic response based upon 
commands; 

starting with a current state of the system, a quadratic programming module formulating a 
problem of achieving the desired dynamic response based upon the current state of the system 
for a window spanning multiple time steps as a solution to a quadratic program problem using 
methods of model predictive control; 

a quadratic programming solver solving the quadratic programming problem in each time 
step using an iterative algorithm in which a Large Sparse Matrix Equation, LSME, is formed 
based upon the model predictive control problem and a quadratic programming algorithm, the 
solver forming a square root of H, Hr, in a first iteration of the iterative algorithm, the solver 
forming a square root of a Large Sparse Matrix of the LSME, LSMroot, using Hr in a quadratic 
programming solver in iterations including the first iteration and a plurality of iterations 
subsequent to the first iteration of the iterative algorithm, the solver completing the solution of 
the LSME based upon LSMroot in iterations including the first iteration and the plurality of 
iterations subsequent to the first iteration of the iterative algorithm , the solver outputting a 
plurality of actuator commands for achieving the desired dynamic response based upon the 
solution to the quadratic programming problem to control a plant. 

19. (Original) The system of claim 18 wherein the solver forms Hr in the first 
iteration and does not recalculate Hr in the plurality of iterations subsequent to the first 
iteration of the iterative algorithm. 

20. (Original) The system of claim 19 wherein the solver forms Hr without first 
forming H. 
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2 1 . (Original) The system of claim 20 wherein H = Hr' * S *Hr and Hr = L W~ where 
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22. (Original) The system of claim 2 1 where the submatrices F, G, and M are 
produced by the recursive sequence of QR factorizations: 
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for various future controller update time points in the prediction horizon. 



23. (Original) The system of claim 18 wherein the solver forms LSMroot without 
first forming the LSM. 
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24. (Original) The system of claim 23 wherein the solver forms the LSMroot by 
solving a normal or generalized QR factorization. 



25. (Original) The system of claim 23 wherein the solver forms the LSMroot by 
solving a generalized QR factorization customized for an interior point method. 

26. (Original) The system of claim 25 where the generalized QR factorization is 



[°] = f ,.«, where «n=[7].uf Ju = [' 



where Pin and S are inputs and triangular P and, possibly, U are outputs, and where 
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27. (Original) The system of claim 23 wherein the solver forms LSMroot by solving 
a generalized QR factorization customized for an active set method. 

28. (Original) The system of Claim 27 where the generalized QR factorization is 

8 V2p = S V2 V . where pin = j-TjT ^ jjT s _jj = s 

where Pin and S are inputs and P and Sr, and possibly, U are outputs and where 



29. (Original) The system of Claim 28 where the generalized QR factorization is 
accomplished via a standard Householder transform based QR algorithm that is modified by 

( vv r l 

replacing the sequence of Householder transforms U i = \ I -2 — — , which are normally 
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used with a corresponding sequence of generalized transformations U i = 




U is unneeded but, if desired, is the cumulative product of the Ui. 

30. (Original) The system of claim 29 where the result of generalized QR 
factorization, which is a triangular matrix consisting of rows containing either purely real 
numbers or purely imaginary numbers, is further factored into the product of a diagonal 
matrix, Sr 2 , which has either 1 or the square root of -1 on the diagonal, and a purely real 
triangular matrix and wherein a square of this diagonal matrix, Sr, will then have either 1 or - 
1 on the diagonal, and Sr and P are purely real numbers. 

3 1 . (Original) The system of claim 27 wherein the solver factors the LSM to use only 
real numbers as: 



where S and SR are constant diagonal matrices with only 1 or -1 entries along the 
diagonal, and where LSMroot is a block triangular matrix where each diagonal 
submatrix is also triangular. 

32. (Original) The system of claim 18 wherein the LSME has the form: 



where z is a vector containing states, controls and state equation adjoint variables for 
each time point, grouped by time; m includes adjoint variables for inequality constraints, 
grouped by time; f and K are vectors; H and J are banded matrices; and T is a diagonal matrix or 
zero. 
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33. (Cancelled) 



34. (Currently Amended) The method of claim 1 further including the step of 
operating the a-plurality of actuators on the a -plant based upon the plurality of actuator 
commands, wherein the feedback sensors are measuring physical variables of the plant. 



35. (Cancelled) 



36. (Currently Amended) The system of claim 1 8 further including a plurality of 
actuators acting upon the plant system based upon the plurality of actuator commands to effect 
the desired dynamic response. 
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